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Vortex oscillations in confined Bose-Einstein condensate interacting with ID optical 
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We study Bose-Einstein condensate of atomic Boson gases trapped in a composite potential of 
a harmonic potential and an optical lattice potential. We found a series of collective excitations 
that induces localized vortex oscillations with a characteristic wavelength. The oscillations 
might be observed experimentally when radial confinement is tight. We present the excitation 
spectra of the vortex oscillation modes and propose a way to experimentally excite the modes. 
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1. Introduction 

Since the realization of vortex state in Bose-Einstein 
condensates (BEC) in traps, 1,2 various kinds of vortex 
configurations are experimentally created. A clear ad- 
vantage of studying vortices in the condensates of atomic 
gases is that the vortex cores are directly observable. Us- 
ing the imaging of atom density in planes perpendicular 
to the vortex axis, vortex lattices, 3 giant vortex, 4 core- 
less vortex, 5 and doubly quantized vortices 6, 7 have been 
observed. By the imaging parallel to the vortex lines the 
bending 8 and the oscillation 9 of the vortex line are found. 

Recently dilute atomic gases interacting with the opti- 
cal lattice 10 are actively studied. The Boson gases in the 
optical lattice behave in two different ways, as a super- 
fluid or a Mott-insulator. 11 The condition of the phase 
transition between these two phases is determined by 
various parameters, such as the particle number and the 
intensity and the dimension of the lattice potential. Vor- 
tices are also an interesting phenomena here. Vortices 
may exist in presence of the optical lattice potentials, 
and the properties of the vortices are affected by the lat- 
tice potential. There also are researches on several kinds 
of the vortices on the BEC in the optical lattices near 
the Mott-insulator transition 12 and in the superfluid re- 
gion. 13 ~ 15 

The vortices whose cores are perpendicular to the 2D 
optical lattices are studied in ref. 13. The dynamics of 
vortices and a new kind of vortex that is unique to the 
optical lattices are reported. On ID optical lattice, the 
vortices parallel to the optical lattice are studied 14, 15 and 
they predicted various vortex oscillations utilizing the an- 
alytical formulations. The formulation assumes periodic 
system along the lattice and neglects the effects of the 
axial confinement potential. 

In this paper, we treat the condensate in ID optical 
lattice confined along the principal axis, to find features 
unique to the confined systems. The whole confined sys- 
tem is treated within the Gross-Pitaevskii (GP) and Bo- 
goliubov equations, and we study structures of the spec- 



tra that are unique to the confined system. 

In our previous work, we found 16 characteristic excita- 
tions in the condensate confined in ID optical lattice in 
absence of the vortex. The branch of these characteristic 
excitations was identified by a breaking down of phase 
count, but the detailed origin and the properties were not 
clarified. This paper studies corresponding modes, espe- 
cially within the Kelvin modes that accompanies vor- 
tices. And identify the characteristic excitations as the 
shortest wavelength modes both in the vortex systems 
and the vortexfree systems. 

2. Formulation 

We study the excitation spectra of BEC confined in a 
composite potential 

V(r, z) = ^{2-KVrfr 2 + V z {z) + V opt {z), (1) 



V z (z) = ^(2nv z fz\ 
V pt(z) = S-E R sin 2 (n^J 



(2) 
(3) 



that consists of an optical lattice in the principal direc- 
tion of consideration and harmonic confinement poten- 
tials along the lateral dimensions. Here s is an inten- 
sity parameter of the periodic optical potential, = 
jjjpj (§) is a recoil energy, m is the mass of a 87 Rb 
atom, v r and v z are the frequencies of the harmonic po- 
tentials, and d is the period of the potential. We em- 
ploy d = 0.395 /urn which is one-half of the typical 17 
wavelength A for the laser beams, v r = 100s _1 and 
v z = 30s _1 which are also within the typical ranges in 
experiments. 

The condensate may have finite angular momentum 
which is essential to treat the systems having vortices. 14 
For simplicity, we assume the condensate rotational sym- 
metry 



(r, z)e l 



(4) 
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Fig. 1. Displacement of the vortex core at the centers (z = jd, 
j is an integer) of the optical lattice sites at s = 5. (a) - (d) 
c = to 3, (e) c = 64, (f) c = 51, and (g) c = 36. 



where w = for vortexfree system and w = 1 for systems 
with a vortex. The wavefunction 0(r, z) is a real function. 
If this function is complex and has vorticity in rz plane, 
the condensate will have parallel vortex rings 18 which is 
also an interesting subject. 

The BEC in an optical lattice with the lower value of 
the intensity s is well described by the Gross-Pitaevskii 
equation 19 

(-CV 2 + V(r, z) + g\4(r)\ 2 ) 0(r) = ^(r) (5) 

where <f>(r) is the wavefunction of the condensate, \i is 
the chemical potential, g = Avh a is an interaction pa- 



rameter, and C — fr? /(2m) is a constant. We denote 
the components of r to be (r,z,6). We treat a conden- 
sate with 5 x 10 4 atoms of 87 Rb. The scattering length 
a = 5.4 nm is employed for the Rb atoms. The period 
of the optical lattice is d = 0.395 ^m and the conden- 
sate occupies about 80 sites. The centered lattice site 
(-d/2 < z < d/2) is occupied by about 1290 (1110) 
atoms when s = (s = 10). 

The condensate supports excitation spectra within the 
Bogoliubov framework. The Bogoliubov equations 19 



(-CV 2 + V(r,z) + 2. 9 |0| 2 - M ) 



(6) 



(-CV 2 + V(r, z) + 2g\0\ 2 - M ) v q - g<P* 2 u q = -e q v q 

(7) 

yield the excitation energies e q and the corresponding 
wavefunctions, u q and v q . The symmetry of system en- 
ables us to rewrite the wavefunctions 



u q (r)=u q (r,z)e^ w+ ^ e , 
v q (r)=u q (r,z)e^ w -^ d , 



(8) 
(9) 



where u q (r, z) and v q (r, z) are real functions and qg is an 
integer representing the angular momentum. The excita- 
tions with angular momentum qg = — 1 is related to the 
oscillation of vortex core and are called Kelvin modes. 



3. Vortex Oscillations 

We consider the small perturbation around the sta- 
tionary state 

<P'(r,t) = ( (>(r) + C 1 (u q (r)e- is ^ h - t£(r)e fe «*/ ft ) (10) 

where C\ is a constant. At t = 0, the perturbation by 
the Kelvin modes (qg = — 1) is 

41 {r, 0) = 0(r, z)e w + & (u q (r, z) - v q (r, z) e - 2lf> ) (11) 

and the corresponding modified density n! = |(/)'| 2 at 
t = is 

n'(r, z, 9) ~ n(r, z) + 2Ci<j>{r, z) 

x {u q (r, z) cos8 — v q (r, z) cos3#} (12) 

where n = \4>\ 2 - The principal direction of perturbation 
of the density n'(r,z,6) above is along the a;-axis (9 = 
0and7r). We write the displacement of the vortex core in 
the xz plane Ax(z), which is defined at bottoms of the 
lattice sites z = jd (j is an integer). Figures 1(a) - (g) 
show several patterns of the displacement Ax(z) of the 
vortex core. 

To balance the original density and the perturbation, 
we employ a value 

max r {n(r)} 



Ci =0.1- 



(13) 



m«jcg, r [\4>(r) {u q (r) - v q (r)} |] ' 

in the figures. Here q covers the Kelvin modes (qg = —1). 

On unconfined periodic system, a wavenumber, e.g. k, 
denotes the phase of the wavefunction through a factor 
exp(zfcz). But explicitly, this definition of k cannot be 
used in confined systems. The exception is the conden- 
sate interacting with moving optical lattices. The veloc- 
ity of the optical lattice and the wavenumber are closely 
related, and the wavenumber remains good approxima- 
tion to denote a excitation. 10, 20,21 Because the moving 
optical lattice is not interest in this article, we employ 
another index similar to the wavenumber, but better re- 
lated to an observable property of the excitation. 

We define a number c that counts the change of sign of 
Ax. Because v q (r, z) = at r — 0, the sign and the ampli- 
tude of Ax corresponds to those of u q (r, z) at r = 0. Fig- 
ure 2(a) plots displacement count c vs. excitation energy 
e for system with (bullets, circles) and without (crosses) 
the optical lattice potential. As the optical lattice poten- 
tial becomes stronger (s increases), the excitation energy 
for certain c becomes lower. It means that vortex oscil- 
lations with many nodes are more likely to be excited at 
stronger lattice potential. 

4. Shortest wavelength modes 

The displacement count c of the Kelvin modes 
smoothly increase as shown in Fig. 2(a). But this in- 
crease breaks down and the c decreases at higher energy, 
see Fig. 2(b). The largest displacement count is around 
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Fig. 2. (a)(b) Displacement count c of the vortex core for s = 
(x), s = 5 (o), and s = 10 (•). (a) The upper horizontal dotted 
line indicates half of the energy of counter-rotating quadrupole 
mode [see Fig. 3(f)]. The dashed horizontal line shows the en- 
ergy —e/2 where e here is the energy of the lowest co-rotating 
quadrupole mode [Fig. 3(a)]. (b) The same plot with wider ranges 
of c and e. The modes whose wavefunction u q has node along 
the radial direction arc omitted. 



c ~ 60, which is close to the number of the occupied 
sites. 

Excitations above the breaking down energy have com- 
mon features in their wavefunctions and the correspond- 
ing oscillation patterns. The oscillation ranges whole the 
condensate [Fig. 1(e)] at the breaking down, but at the 
higher energies the z ranges of the oscillations shrink. 
Accompanying this shrink, the displacement count c de- 
crease to, for example, 36 [Fig. 1(f) - (g)]. 

These plots show that the displacement Ax in this de- 
creasing aspect has opposite sign between the neighbor- 
ing sites. It means that the wavefunctions u at the vortex 
core (r = 0) have opposite phase between the neighbor- 
ing sites. Therefore, the wavelength of these excitations 
are 2d. This is the shortest wavelength possible within 
the range of excitation energy e <C sEr we consider. We 
call these excitations the shortest wavelength excitations 
hereafter. 
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Fig. 3. Wavefunctions u(r, z) of the quadrupole (qg = +2, —2) 
modes at intensity s = 5. Solid (dotted) lines show contours at 
u(r,z) > 0.01 (u(r,z) < —0.01). Dashed lines are Thomas-Fcrmi 
profile of the vortexfree condensate without the optical lattice, 
(a) - (b) First two excitations within qg = +2. (c) - (d) Two 
excitations with highest energy within the first band (qg = +2, 
radial node free) .Note that each of these two have opposite sign 
(i.e. opposite phase) of u between the neighboring sites, (e) The 
first excitation within qg = —2. (f) The first centered excitation 
of qg = —2. Vortexfree [w = in eq. (4)] condensate supports 
quadruple (qg = ±2) modes similar to these (a) - (d). 



We pointed out the shortest wavelength excitations 
within the Kelvin modes. In addition to these, excita- 
tions of this kind commonly exist in the confined con- 
densates interacting with an optical lattice. The conden- 
sate with and without the vortex supports shortest wave- 
length excitations within the breathing (qg = 0), dipolc 
(qg = ±1), quadrupole (qg = ±2) modes, and modes with 
several higher angular momenta. Figures 3(c) and 3(d) 
are those within the quadrupole modes of qg = +2. One- 
dimensional system with ID optical lattice and harmonic 
confinement potential also has corresponding excitations. 
The set of these excitations are mentioned as "branch B" 
in ref. 16. 

4-1 Energy of the shortest wavelength modes 

In an unconfincd 1-dimensional periodic system with 
period d, the largest wavenumber is ir/d, and the corre- 
sponding phase difference between the neighboring sites 
is 7r at the upper end of the so-called first band. The 
shortest wavelength modes in our confined systems also 
have the phase differences n. And there is always an en- 
ergy gap above the shortest wavelength excitations if we 
ignore the modes with radial nodes. So it is reasonable 
to think that these excitations also belong to the upper 
end of the first band. 

Figure 4(a) compares ranges of the energy of the short- 
est wavelength excitations and the upper end Em(s) of 
the first band of an ideal (noninteracting, unconfincd) 
1-dimensional periodic system. 10 ' 16 ' 21 It is seen that 
the energy ranges of the shortest wavelength excitations 
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Fig. 4. (a) Boxes shows ranges of energy (vertical) of the Kelvin 
modes and s e g (horizontal) for systems with lattice intensities 
s = 5,7, and 10. Solid curve is energy E±d (s) of upper end of the 
first band for the nonintcracting ID periodic system, (b) Ener- 
gies i?iD( s eff( 2 )) corresponding to the local value of the effective 
lattice intensity s e ff (z) and spatial distributions of wavefunctions 
of the Kelvin modes are compared. Solid curve is i?iD(s e ff( z ))- 
Sets of the horizontal short lines are the wavefunctions u q (r,z) 
at r = at its energy e q . Bold (thin) lines shows u q (0, z) is pos- 
itive (negative). Axial (z-axis) distributions of the Kelvin modes 
are localized near the peak of Ei-£,(s e g(z)). Both of the arrows 
in (a) and (b) shows the same range of -EiD^eflK 2 )) f° r s = 5. 



overlap with Em(s). The energy ranges of the short- 
est wavelength excitations of another angular momenta 
qe = 0,1, ±2 also overlap Em{s) despite the change in 
the angular momenta. So the difference of angular mo- 
mentum is not essential in the comparison in Fig. 4(a). 

Range of Em is < Em < £r (numerically calculated 
from eq. (12) in ref. 16). So the ranges of the shortest 
wavelength modes for various s will also be within < 
Em < En. 

4-2 Origin of the spatial dependence 

While other modes are elongated whole the conden- 
sate, he shortest wavelength excitations are localized 
near the center (z = 0) of the system. This feature is 
clearer for modes with the highest energies, e.g. Figs. 1(g) 
and 3(d). 
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Fig. 5. Displacement count c vs. energy of the Kelvin modes for 
s = 5 (x) and s = 9 (•). The trapping potential parameters are 
v r = 4000 Hz, v z = 300 Hz. The horizontal lines are the same as 
Fig. 2(a). 



Within the GP equation and the Bogoliubov equa- 
tions, the depth of the optical lattice s is changed due to 
the meanfield energy g\(j)\ 2 in eqs. (5), (6), and (7). The 
effective depth becomes smallest near the center (z = 0) 
of the system. In ID periodic system, the smaller s be- 
comes, the higher the upper edge Em(s) of the first band 
becomes. This agrees with the localization of the high- 
est shortest wavelength mode around z = at which the 
effective s becomes the smallest. 

To make this discussion clear, we define effective am- 
plitude 



Sctt(z) =s- max (g\<j>(r, z)\ 2 ) /E R 



(14) 



at bottoms of the lattice sites z = jd (j is an integer). 
Within the 'max' above, the radius r is chosen to have 
the maximum. The function s e ff(z) is smallest near the 
center (z = 0) of the system. 

Spatial (along z-axis) dependence of the effective 
depth s c ff(z) and the energy may be related by a func- 
tion Em(s e s(z)). This function and the energy and spa- 
tial localization of the Kelvin modes are compared in 
Fig. 4(b). The excitations at the left hand side above 
e > 0.27 En are the wavefunctions of the shortest wave- 
length modes. 22 The correspondence between the s e ff and 
the axial (z-) range of the wavefunctions is clear. 

5. How to observe the shortest wavelength 
modes 

On trapped BEC of atomic gases, the quadrupolc 
modes may be excited experimentally 9 by deforming the 
confining potential. And both the experiment 9 and the 
following theoretical study 23 insist that Kelvin modes 
(qg = —1) are excited by the decay from the counter- 
rotating quadrupole mode {qe = —2). For the conden- 
sates interacting with the periodic potentials, a reso- 
nance condition between the quadrupole mode and the 
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Kelvin modes is given in § V of ref. 14. The criteria given 
in the references above may be simplified to this: The ex- 
citation energy e_i (or frequency E-i/h) of the excited 
Kelvin modes satisfy 



£-1 = £-2/2, 



(15) 



where e_2 is the energy of the counter-rotating 
quadrupole mode at qg = —2. If one of the shortest wave- 
length modes satisfies eq. (15), the mode will be excited 
as an oscillation of the vortex core. The oscillation will 
be experimentally observable. 

The corresponding excitations in this paper is the low- 
est centered counter- rotating quadrupole mode qg = — 2 
in Fig. 3(f) and the Kelvin mode whose energy satis- 
fies eq. (15). The excitation energy £_ 2 varies only little 
(stays between 0.0343£ R to 0.0351£r) for s = to 10. 
And half of the energy £-2/2 is indicated in Fig. 2(a). 
The figure shows that Kelvin modes with c = 13 (c = 22) 
at s = 5 (s = 22) are on the dotted line and satisfying 
the criterion. 

Then how can we make the shortest wavelength modes 
satisfy eq. (15)? The right hand side, the energy of the 
centered counter-rotating quadrupole modes is e_ 2 — 
1.2 hv r where h is the Plank constant. The left hand 
side, the energy of the shortest wavelength modes, ranges 
< e < -Er from the discussion in § 4.1. As _Er, is a 
function of period d of the optical lattice and £-2 is ap- 
proximately a function of radial confinement potential 
parameter v r , these two may be varied independently 

We increased v r while keeping d. At v r = 4000 Hz, 
the shortest wavelength mode satisfies eq. (15). Here ax- 
ial trap frequency v z particle number N are changed to 
v z = 300 Hz and N = 10 5 . Figure 5 presents spectra of 
the Kelvin modes and the energy £-2/2. In this tight ra- 
dial trapping potential, the shortest wavelength modes 
satisfies eq. (15) at s > 5 and localized vortex oscillation 
of period 2d will be observed by exciting the counter- 
rotating quadrupole mode. 

In writing eq. (15), we neglected the width of the res- 
onance whose evaluation 23 requires huge computational 
resource. Explicitly speaking, the condition is satisfied at 
only discrete values of s. The allowed width of s is a func- 
tion of various parameters, including the temperature. 
This discreteness can be avoided by smoothly varying s. 

6. Discussion 

In BEC of atom gases in ID optical lattice, existence 
of the spatially localized modes was reported in ref. 16. 
But the nature of the modes was not clarified there. In 
this paper, we find that the modes have common wave- 
length 2d and commonly exists at various angular mo- 
menta, and exists both at the vortexfree condensate and 
the condensate with a vortex. 

Among the modes, we report on the Kelvin modes and 



the corresponding vortex oscillations in systems with vor- 
tex. Provided that the vortex oscillation was observed 9 
in systems without the optical lattice, the same obser- 
vation will be possible with optical lattice. By choosing 
the trapping potentials and the intensity of the optical 
lattice, we can make the vortex oscillation modes of wave- 
length 2d satisfy the resonance condition eq. (15). Then 
the corresponding oscillation, whose directions are the 
opposite between the neighboring sites, will be observed. 

I thank S. K. Yip and T. Mizushima for useful discus- 
sions. 
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